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We present a toy model that exhibits clash-of-symmetries style Higgs field kink configurations 
in a Randall-Sundrum-like spacetime. The model has two complex scalar fields $1,2, with a sextic 
potential obeying global U(l)(g)U(l) and discrete $1 <-> $2 interchange symmetries. The scalar 
fields are coupled to 4+ 1 dimensional gravity endowed with a bulk cosmological constant. We show 

in ■ 

that the coupled Einstein-Higgs field equations have an interesting analytic solution provided the 
sextic potential adopts a particular form. The 4+1 metric is shown to be that of a smoothed-out 
' Randall-Sundrum type of spacetime. The thin-brane Randall- Sundrum limit, whereby the Higgs 

^jjy field kinks become step functions, is carefully defined in terms of the fundamental parameters in the 

r—{ \ action. The "clash of symmetries" feature, defined in previous papers, is manifested here through 

the fact that both of the U(l) symmetries are spontaneously broken at all non-asymptotic points 
in the extra dimension w. One of the U(l)'s is asymptotically restored as w —* —00, with the 
other U(l) restored as w — » +00. The spontaneously broken discrete symmetry ensures topological 
stability. In the gauged version of this model we find new flat-space solutions, but in the warped 
' metric case we have been unable to find any solutions with nonzero gauge fields. 
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The breaking of symmetries is one of the most important issues in particle physics. Nature has taught us the 
relevance of symmetry principles, especially through Noether's connection of global symmetries with conservation 
laws and the successful description of fundamental interactions via local gauge invariance. But to reconcile symmetric 
field equations with the real world, one needs that some of the solutions to those equations have a reduced symmetry 
compared to the equations themselves. This is, of course, what we mean by 'spontaneous symmetry breaking (SSB)' 
. or 'hidden symmetry'. 

The precise mechanism by which electroweak SSB happens has yet to be experimentally determined. In the standard 
model (SM), the default option is to use an elementary scalar or Higgs field doublet with a quartic potential that 
induces a nonzero vacuum expectation value for one of the components. This will soon be tested at the Large Hadron 
Collider through its physical- Higgs boson search programme. 

The purpose of this paper is to further develop a different method of symmetry breaking called the 'clash of 
symmetries' Q, 0, El Q (see also Q). This mechanism is relevant for brane- world extensions to the regular 3 + 1 
dimensional standard model 0, Q • The important development we report here is the rigorous incorporation of gravity 
into one of the toy Higgs models explored in the previous papers on this topic [3j . Since gravity is fundamental to 
brane-world models, this is a very welcome step forward. 
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The toy model we shall analyse consists of two complex scalar helds $1.2 coupled to 4 + 1 dimensional gravity 
endowed with a bulk cosmological constant. The action is constructed to be invariant under global U(l)(g>U(l) 
transformations and the discrete <&i <-> $2 interchange operation. The Higgs potential we construct is sextic, and 
its specific form is driven by our desire to find an analytic solution to the coupled Einstein- Higgs field equations. 
The spacetime metric tensor we posit has the Randall-Sundrum non-factorisable form but with the warp factor 
a completely smooth function of the extra dimension coordinate. The precise Randall-Sundrum metric emerges by 
taking a certain limit. The Higgs field configurations we seek are kink-like, associated with the spontaneously broken 
discrete symmetry. 

We also study the gauged version of this model, extending the analysis of Ref.|3j]. We find additional non-trivial 
flat-space solutions, but despite an extensive numerical search we have been unable to find any warped metric solutions 
with nonzero gauge fields. 

Section ITT1 reviews the clash of symmetries mechanism, while Sec. IIIII defines and analyses the ungauged version of 
our toy model. This is followed by an account of the gauged cases in Sec. IIVI The final section is a conclusion and 
outlook for future work. 



II. THE CLASH OF SYMMETRIES 

The default SSB mechanism employed in the SM and many extensions to it utilises a simple spatially homogeneous, 
static and stable solution to the field equations: the Higgs fields <j)i(x^) are set equal to those constant values 
corresponding to an absolute minimum of the Higgs potential. The locations of the nonzero constants within the 
Higgs multiplets then determine the stability group of the vacuum state, what we commonly call the unbroken 
symmetry after SSB. The nonzero Higgs fields become constant background scalar fields throughout the universe. 

There are other types of solutions to Higgs field equations that can also serve as stable, static background fields: 
topological solitons, such as kinks, strings and monopoles. While there is much interest in such configurations, they 
are not usually considered as viable candidates for playing a major role in spontaneous symmetry breaking. This is 
simply because they are spatially inhomogeneous. The energy densities required for electroweak and higher symmetry 
breaking are so high as to be incompatible with the strong evidence for a universe that is spatially homogeneous at large 
scales. This objection, however, does not apply to brane-world models, because the non-trivial spatial dependence 
can be restricted to the extra dimension coordinates only. 

For definiteness, consider the case of one extra dimension described by Cartesian coordinate w and topologically- 
stablc Higgs configurations <j>i(w), some of which have kink-form with respect to w. The pattern of spontaneous 
breaking then effectively becomes a function of w. The 3 + 1-brane is located at, say, w = 0. If the degrees of freedom 
that are confined to the brane [Eo| are confined absolutely, then the unbroken symmetry is the stability group of 
4>i(w — 0), to be denoted by H(w = 0) from now on (and let the internal symmetry group of the action be G). The 
symmetry breaking will be communicated to the brane-world denizens through appropriate interactions. But in a 
quantal (and perhaps even in a classical) world, one would not expect strict confinement to the brane: there should 
be some leakage off it into the bulk. One therefore expects the brane-localised fields to also couple, but with reduced 
strength, to <j>i(w) states where w is slightly different from zero. If the stability group of 0j(O < |w| < e), call it 
H(w ~ e), is different from H(w — 0), then a rather rich effective symmetry breaking outcome could ensue for the 
brane world. For example, a hierarchical breaking pattern G — > H(w = 0) — > H{w ~ c) was speculated upon in 
Ref.Q, where in that case G = SO(10), H(w = 0) = SU(4)®SU(2)®U(1) and H(w ~ e) = SU(3)®SU(2)®U(1) 2 . 
The potential application to grand unified models and the gauge hierarchy puzzle is clear. 

The 'clash of symmetries' phenomenon allows the unbroken symmetry at non-asymptotic points in the extra di- 
mension, I to I < 00, to be smaller than the symmetry holding at \w\ =00 0, E3> Q] ■ The boundary conditions at 
infinity for kink configurations are set to be the Higgs vacuum expectation values (VEVs). They force the spontaneous 
breaking G — > i? — * 00). In usual SSB, the Higgs VEV configuration is used throughout the whole of space, so the 
symmetry is broken to the group we are calling -ff (|w| — > 00) everywhere, homogeneously. In the clash of symmetries 
mechanism, this symmetry breaking pattern only holds asymptotically with respect to the extra dimcnsion(s). 
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The clash of symmetries can occur if the isomorphic subgroups left unbroken at w — — oo and w = +00, H{w — > 
—00) = if (— 00) and H(w — > +00) = if (+00) respectively, can be differently embedded within the parent group G. In 
the cases examined so far in the literature, the breakdown at finite values of w is to the intersection of the asymptotic 
stability groups, 

H(\w\ < 00) = H(-oo) n H(+oo) = H clash . (2.1) 

Although H(—oo) and if (+00) are isomorphic, their different embeddings within G leads their intersection to be a 
smaller group. 1 Some special values of w, most often w = 0, may correspond to stability groups larger than if c i a sh 
because the spatially varying Higgs configurations may instantaneously pass through special patterns as w varies. 
For instance, a component of a Higgs multiplet may be an odd function of w such as tanh(w), so it would vanish at 
w = and thus the unbroken symmetry would be enhanced at that point. This possibility lies behind our remarks 
two paragraphs above about an application to the gauge hierarchy problem. 

Discrete symmetries also play a vital role: their spontaneous breakdown leads the vacuum manifold to have discon- 
nected pieces. Kinks interpolate between vacuum states from disconnected pieces. Within a given topological class 
of kinks, the one with lowest energy is then guaranteed to be (topologically) stable. 

The easiest example of the clash of symmetries uses G = SU(3) spontaneously breaking to SU(2). An SU(2) 
subgroup can be embedded in three different ways in SU(3): f-spin, [/-spin and V-spin. Clash of symmetries kinks 
interpolate between VEVs respecting (i, U), (U, V) or (V, I) at w = (—00, +00). The differently embedded asymptotic 
SU(2) symmetries "clash" at non-asymptotic points, leading if c iash to be SU(2)/nSU(2)[/ = {1} and so on. 2 

The model-building promise of the clash of symmetries lies in the greater symmetry breaking power of spatially 
inhomogeneous Higgs configurations over the simple homogeneous alternative of standard SSB. We hope that this will 
lead to models with simpler Higgs sectors and fewer associated parameters. Some of us have already speculated that 
the ultimate application might be to an Eg model, because of the fact that there are three different Eq generators that 
can serve as electric charge, corresponding to three different embeddings of electromagnetism within the full group 
0. Standard Eq models can only use one at a time. Even more speculatively, a connection between that threefold 
structure and the three families of quarks and leptons might exist 3 Q]. 

But many preparatory studies are needed before such an ambitious vision can be realised. In this paper, we take 
another important step along that road: the incorporation of gravity. To facilitate the study of gravity, we begin by 
simplifying the group theoretic structure as much as possible without losing the essence of the clash of symmetries 
mechanism: we take G = U(l)i(g>U(l)2, with if(— 00) = U(l)i and ff(+oo) = U(l)2 leading to ff c iash = {1}< The 
discrete symmetry we need for topological stability interchanges the two Abelian sectors. 

In the next section we show that the ungauged version of this model can yield a clash-of-symmetries-style kink in 
a warped metric spacetime. The section after that revisits the gauged extension in both flat and curved spacetime. 



1 There will also typically be kink configurations that interpolate between identically embedded subgroups. Such kinks do not display the 
clash of symmetries, and one must ensure that they have higher energy than those displaying the clash. 

2 To actually implement this simple SU(3) scheme requires more work: see Ref. for details. 

3 There is some similarity between the clash of symmetries mechanism and orbifold symmetry breaking 0. In particular, Ref. 
constructs a model where different subgroups of SO(10) are left unbroken on different branes, with the overall unrboken symmetry being 
the intersection of those individual subgroups, and with families also distributed amongst the branes. While there are similarities, it 
is clear that the mechanisms are not identical, because the clash of symmetries has a physical Higgs field while the orbifold method is 
Higgsless. 
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III. THE UNGAUGED MODEL 



A. The field equations 



The model has two complex Higgs fields $! and $2 coupling to 4 + 1 dimensional gravity, and it respects the 
$1 <-> $2 discrete symmetry. The action is 



S=- 



R + g 



MN 



(dM<S , id N <Z> 1 +d M $ 2 d N <i>2) + V($i$i,$ 2 $2) V^9 db ^ 



(3.1) 



where <?mjv is the 4+1 dimensional metric tensor, with M, N = 0, 1, 2, 3, 5, R is the curvature scalar, and overbars 
denote complex conjugation. The matter action S m is defined as S with the Einstein-Hilbert contribution absent. 
The Higgs potential V is as yet unspecified except that it must be invariant under $1 <-> $2- It will have a constant 
piece that serves as the bulk cosmological constant. We employ the (— , +, +, +, +) signature and our sign conventions 
are as in Weinberg (apart from our definition of g = Det(<?)) |ll| . 

The Einstein equation found by requiring a stationary action under metric variations is 

1 



G MN = _^ T 
K 



MN 



(3.2) 



where the Einstein tensor is, as usual, G MN 
tensor, 



R MN -(l/2)g MN R with R MN being the Ricci tensor. The stress-energy 
2 SS m 



T 



MN 



-9 Sg M N 



evaluates to 



where 



^N = 2g MP g NQ tpQ _ g MN {g PQ tpQ 



t MN EE d M $i d N $! + d M ®2 d N § 2 - 



The Higgs Euler-Lagrange equations are 







(3.3) 

(3.4) 
(3.5) 

(3.6) 



where i = 1,2, and the complex conjugate equations also hold. 

We now specify our solution ansatz. We seek configurations that only depend on the extra dimension x 5 = w. The 
imaginary parts of the Higgs fields are set to zero and we write 



4>, 



V2 ' 



(3.7) 



where 4n{w) is real. The metric ansatz, of Randall- Sundrum type, is defined through the 4+1 dimensional line 
element 



ds\ = dw 2 + e 2 l {w Usl 



(3.8) 



where f(w) is an as yet unknown function and ds\ = r] fl ^dx ti dx 1 ' is the 3 + 1 dimensional Minkowski line element (our 
brane will be Minkowski), with ^t, v = 0, 1, 2, 3 or t, x, y, z. The field equations then reduce to 



3.r+6(f) 2 = — 

K 



6(/') 2 
- 4/ V 



1 

K 

dV 



\ W? + (0 2 ') 2 ] - V 



(3.9) 
(3.10) 
(3.11) 
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where prime denotes differentiation with respect to w. Adding and subtracting Eqs. I|3.9[l and (|3.10f) gives 

/" + 4(/') 2 = ~V, (3.12) 

/" = ~[(<PiT + (fc') 2 ]- (3-13) 
3k 

Note that Eq. (|3.13|l is independent of the potential, a feature we will use below. 
We will construct a potential V that has degenerate global minima at 

vacuum 1 : <j>\ = u, 4>2 — and vacuum 2 : <p\ = 0, 4>2 — u, (3-14) 

where u would be the VEV in the standard homogeneous SSB mechanism. The vacua are related by the now 
spontaneously broken discrete interchange symmetry. The Higgs field boundary conditions will be 

<j>i(— oo) = vacuum 2 and 0i(+oo) = vacuum 1. (3.15) 



B. Constructing the solution 

In 3 + 1 dimensions without gravity, renormalisability limits us to quartic potentials. With extra dimensions, and 
with gravity included, this is no longer a well-motivated requirement. In fact, the only general constraint on the form 
of V we will impose is that it be bounded from below. A very important topic for the future of the clash of symmetries 
mechanism will indeed be the rules by which potentials should be specified. 

Our more modest goal here is to provide an existence proof that kinks featuring the clash of symmetries form can be 
consistently coupled to Randall-Sundrum-like gravity. We will search for an analytically tractable model for reasons 
of both simplicity and elegance. To do this, it makes sense to invert the problem by first specifying the solutions 4>i (w) 
and f(w) and then constructing a potential that gives those solutions. This, in itself, is non-trivial and interesting. 

It turns out that a good choice, satisfying the boundary conditions of Eq. (|3.15l) and also the Higgs-potential- 
independent Eq. (|3.13|l . is 

u 



<j>l{w) = -^V 1 + tanh/3w, (3.16) 
v2 



(j) s 2 (w) = —j=yjl~ tanh/3w, (3.17) 
V2 

2 

U 

/» = -—In (cosh /? W ), (3-18) 

where (3 is a free parameter that specifies the brane thickness, and the superscript s denotes "solution" . The Higgs 
configurations are the square roots of the most naive choices, and f s (w) is essentially a regularisation of |io|. The latter 
is, of course, the corresponding function in the Randall-Sundrum scenario. The smoothing out of the Randall-Sundrum 
cusp goes hand-in-hand with the existence of Higgs field kinks coupled to gravity, and the brane is dynamically 
generated rather than put in by hand. Note also that 

#H 2 +#(u>) 2 =u 2 , (3.19) 

a feature that will turn out to be very useful. 

The first inkling for what V should be comes from Eq. I|3.12|l . It suggests that 

where U is a function that vanishes at the level of the solution, that is U((f>l(w) 2 , 4>2( w ) 2 ) — 0- Notice that a negative 
bulk cosmological constant, 

A 5 - (3.21) 
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follows from the ansatz: the bulk has anti-de Sitter qualities, as usual in Randall-Sundrum-like models, though in our 
case it is of course not precisely anti-de Sitter. 

More information about U is supplied by the Higgs field equation (|3.11|) . To solve it, one must specify that 



3 k 



where the new function W obeys 



W(0) 



dW 



n.2 



= 0, 
= 0, 



(3.22) 

(3.23) 
(3.24) 



so as not to affect the satisfaction of the field equations, Eq. (|3.12|) and Eq. (|3.13|l respectively. As far as they are 
concerned, W can vanish. The resulting sextic potential, however, would not be bounded from below because the 
coefficient of the first term in U is negative. 

We therefore construct W so that, as well as obeying Eqs. I|3.23[) and l|3.24|l . it leads to a potential that is both 
bounded from below and has global minima at Eq. I|3.14|l . The simplest suitable function is 



W = ( 



Au 2 



• ^) +<>: 



,2 \2 



(3.25) 



where rj and £ are dimensionless parameters. 

The full Higgs potential for the model is specified by Eqs. (|3.20|l . I|3.22|l and (|3.25|) . We will present a global 
minimisation analysis of it in the next subsection. (To reinstate the complex nature of the Higgs fields, the substitution 
fa 2 — ► 2<£>,;<£>j should be made, though for the forthcoming discussion we need only work with the fa.) 

Before doing so, let us identify the Randall- Sundrum limit of the model. We need to identify a limit in which the warp 
factor, e 2 ^", becomes e~ a \ w \ where a is some positive constant. Following Ref . [l^ . we note that (cosh/Su;) - ^ 1 /' 3 ) — > 
e - '™' as (3 — > oo. Writing 



(cosh (3w) 



we see that the Randall-Sundrum warp factor, e -- ^'™', is obtained in the limit 

/3 — > oo, u — * such that u 2 (3 — > finite constant. 



(3.26) 



(3.27) 



We also require r\ and £ to be unchanged when taking the limit. Note that the Higgs field configurations, Eqs. 
(|3.1ti|l and l|3.17|l . become step functions as (3 — > oo, with their magnitude tending to zero because u — > 0. The 
Randall-Sundrum fine-tuning condition between the brane tension and the bulk cosmological constant is a feature of 
our solution. Defining the brane tension via 



Ai, = lim 

RS 



(3.28) 



it is easy to show that 



At = --kA, 



(3.29) 



which is the Randall-Sundrum condition. 



C. Minimisation analysis 



We now identify the parameter space region where our potential has global minima specified by Eq. (|3.14|) . 
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To establish the parameter region where the potential is bounded from below, we look at the sixth-order terms in 
isolation. They are 



y 6 -A(0 2 + ^)[ C(0 2 + 2 )2 _ 4(/) 2 : 



(3.30) 



where A = (3/3 2 /8u 4 )[l + 2w 2 /(9k)] > 0. The substitutions fa = xcos9 and fa = xsin0 reveal that V 6 = Ax 6 (C 
4 sin 2 9 cos 2 9) , from which the condition 



C>1 



(3.31) 



is trivially seen to be the only requirement to produce boundcdncss from below. 

The minimisation analysis is a four-parameter problem. For some of the local minima, both their (fa, fa) locations 
and the values of the potential at those points are complicated functions of the parameters. But there is a smart way to 
evade these unenlightening complications while yielding physically meaningful results: approach the Randall-Sundrum 
limit defined by Eq. i|3.27[) . Our task then becomes a two parameter problem in (77, that is readily visualised. 

To see this, note that in Eqs. I|3.2U[1 . (|3.22|) and (|3.25|l . the combinations 



1 + 



3 k 



and 



U 

3 k 



(3.32) 



frequently appear. In the Randall-Sundrum limit, which includes u — > 0, these combinations simplify to 1 and 3/2 
respectively. In the exact Randall-Sundrum limit, the meaning to be ascribed to the potential is obscure, so we cannot 
work exactly in that limit. Instead, we will approach the limit without reaching it, and write an approximation to 
the potential justified by the above remarks. It is 



V ~V F 



approx 



f3 2 u 4 


£_ 


24k 


2u 2 


3P 2 u 2 


" u 2 
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9k 



3f_ 

2u 4 



2 2 
■P1P2 



Jl „2 / J2 



,2 \2 



(3.33) 



where 



Pi = — 



(3.34) 



is a dimensionless Higgs field. The fact that (3 contributes only to an overall scale, and that k appears solely in the 
constant term, reduces the analysis to a two-parameter problem in (77, £). While k contributes to the potential values, 
it does so equally for all local minima in this approximation, which means that only 77 and £ need be considered when 
comparing the local minima to see which are also the global minima. 

The potential has three varieties, denoted (I, II, III), of possible local minima (or saddle points), 



(I) 


Pi 


= 1, P2 = 


(II) 


pi 


= p1^o 


(III) 


pi 


= P2=0 



pi = 0, P2 = 1 (discrete symmetry breaking), 



(no SSB). 



(3.35) 
(3.36) 
(3.37) 



Which of these are actual local minima depends on (ij, (). We want I to describe the global minima. 

Of the three, it is // that has the complicated parameter dependence. In fact it is quite easy to show, even for the 
full potential of Eqs. (|3~2U|) . $3~Zfy and (|3~25|) . that V(I) < V(III) requires that 



77 > 1 



(3.38) 



which is thus a necessary but not sufficient condition to ensure that I gives the global minima. Sufficiency requires 
an examination of the parameter dependence of V(IT) also. The results are shown in Fig. where the shaded area 
is the Higgs parameter region where / is the global minimum. 
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FIG. 1: (Colour online.) The regions in the (rj,Q plane where the approximate potential has the required discrete symmetry 
breaking global minimum. The potential is bounded from below for f > 1. In the red (or medium-shaded) region, 4>\ — 02 = 
is a local maximum, <j>\ — <j)% 7^ is a saddle point, and <j>\ — it 2 , 02 = and 0i = 0, 0§ = u 2 are degenerate global 
minima. The green (or light-shaded) region differs from the red (medium-shaded) only in that 0i — 02 = becomes a local 
minimum. As r\ decreases further, the 0i = 02 = minimum keeps falling and in the blue (or dark-shaded) region it lies 



below the 



7^ saddle point. For r) < 1, 0i 



becomes the global minimum. In the unshaded region in the 



77, £ > 1 quadrant, <j>\ = 02 7^ is the global minimum. The border between the shaded and unshaded regions is a smooth but 
complicated function ("(77) given in the main text. 



In the red (or medium-shaded if viewing in black-and-white) region, / is a local minimum and II a saddle point 
with V(I) < V(II), while 777 is a local maximum. As 77 decreases below 3/2, III develops into a local minimum 
also; the ordering is V(I) < V(II) < V(III) in the green (light-shaded) region and V(I) < V(III) < V(II) in the 
blue (dark-shaded). For 77 < 1, III is the global minimum so there is no SSB. The border line between the shaded 
and unshaded areas is given by 

-27 - 18 <q + 61 rf + y/729 + 972 77 - 2970 if + 1900 r? 3 - 375 ^ 



C 



24 if 



(3.39) 



Lowering £ has the effect of lowering V(II). Below this border line, V(II) < V(I) < V(III), so the symmetry 
preserving vacuum is the global minimum. 

The border between the blue (dark-shaded) and green (light-shaded) regions is given by 



V 



8C-3C 2 + 2V3y/4C 2 -7C 3 + 3C 4 
3C 2 



(3.40) 



Figures T2I4I show contour plots of the approximate potential for representative points in the red (medium-shaded), 
green (light-shaded) and blue (dark-shaded) areas of Fig. ^ respectively. The colour- or shade-coding is explained in 
the captions. The plots illustrate the pattern of extrema described above. 



IV. THE GAUGED MODEL 



It is natural to extend our toy model by gauging U(l)g)U(l), and to study the dependence of the gauge fields on 
w by self-consistently solving all of the coupled Euler-Lagrange equations. In the flat space case, Rozowsky, Volkas 
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FIG. 2: (Colour online.) Contour plot of the approximate potential for r\ — £ = 1.6 to illustrate the qualitative features of the 
extrema for a representative point in the red (medium-shaded) region of Fig. Q Within this plot, red colour codes for minima 
and blue for maxima. (These look similar in black-and-white, but the description below clarifies which is which.) The origin, 
pi = p2 = 0, is a local maximum. The degenerate global minima are at (pi,f>2) = (0, ±1) and (±1, 0), while extrema //, where 
pi = P2, are saddle points (light green or very- lightly shaded). 




-1 -0.5 0.5 1 



FIG. 3: (Colour online.) Contour plot of the approximate potential for r\ = 1.3 and C, = 1.6 to illustrate the qualitative features 
of the extrema for a representative point in the green (light-shaded) region of Fig. Q Within this plot, red colour codes for 
minima and blue for maxima. (In black-and-white, the red region appears darker than the blue in this plot; the description 
below clarifies which is which.) The origin, p\ — p2 = 0, is now a local minimum, but the degenerate global minima are still at 
(pijf>2) = (0, ±1) and (±1,0). The relevant p\ = p\ extrema are again saddle points. 
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FIG. 4: (Colour online.) Contour plot of the approximate potential for r\ — 1.05 and £ = 1.6 to illustrate the qualitative 
features of the extrema for a representative point in the blue (dark-shaded) region of Fig.Q Within this plot, red colour codes 
for minima and blue for maxima. (In black-and-white they appear similar, but the description below clarifies which is which.) 
The origin, pi = p2 = 0, is now a deeper local minimum (yellow, or medium-shaded), having dropped below the p\ = p\ saddle 
points (green, or medium-shaded valleys). The local maxima are at approximately (|pi],|f>2|) = (0.5,0.5). The degenerate 
global minima are at (pi,p2) = (0, ±1) and (±1,0). 



and Wali 3] found some interesting solutions where the gauge field coupled to $i tended towards a linear function 
of w on the side of the wall where $; — > 0, and was exponentially suppressed on the other side (Meissner effect). 
The linearly rising gauge field implies an asymptotically constant magnetic field, so the physical picture was of two 
infinite, uniform sheets (one for each sector) of supercurrent density flowing along the domain wall. In this section 
we will first generalise the Rozowsky et al. flat-space analysis and then consider the curved space problem. 



A. Flat space 



The Lagrangian is given by 
1 



C = --F( vn "F 1MN - -Ff»F 2M N - (I9 M $ 1 )*( J D M $ 1 ) - p M $ 2 )*p M $2) - V, 



where 



D M ®i = {9m ~ ieAiM - ieA 2M )$i, 
D M $2 = {9m - ieAxM - ieA 2M )^2, 

and we take V to be quartic 

V = Ai(<j>i$i + $ 2 ( I > 2 - U 2 ) 2 + A 2 $l ( l > l$2 ( i > 2- 

In the Ai ; 2 > parameter space region, the global minima of V are 

Vacuum 1 : ¥i$i = u 2 , <& 2 = and Vacuum 2 : $i = 0, $2^2 = u 2 . 



(4.1) 



(4.2) 



(4.3) 



(4.4) 
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The Lagrangian is invariant under the discrete symmetry 

$1 <-> $2, A\ <-> A 2 . (4.5) 



This model generalises that considered in Ref.|3( by having e^O. 
The Euler-Lagrange equations are 

D M D M $i = 2Ai +¥j ■* i -u 2 ) + A 2 $ i * i $ J -, 
d M F t MN = -2Im(e$ i D N <f> i + e~$ j D N <f> j ), (4.6) 

where = (1,2) or (2,1). We now specialise to configurations that depend only on w, impose Lorentz gauge 

8m Af 1 = 0, and write $j = R4 e l6i . The equations become, 

R'l = R t [(eA? + eA%)(eA i „ + eA j „)]+2\ l R i {R 2 +R?-u 2 ) + \ 2 R i R 2 , (4.7) 
Af = 2eR 2 (eAt^~eA^) + 2~eR 2 AeA^ + eAt), (4.8) 



e\ = -(eAl + eA^). (4.9) 

Observe from the last of these equations that the A = 5 gauge field components are pure gauge, so they and the phase 
fields Qi decouple from the problem. 

We choose to look for solutions that do not break rotational invariance in the x, y and z directions. This means 
that we can always rotate the x, y, z coordinate system so that the only nonzero component of (Af , A\ , Af) is in, say, 
the a;-direction (this is tantamount to the requirement that the ratios of these components are w- independent). We 
also choose that the same component is nonzero for both i = 1 and i = 2. Equations l|4.7|l and l|4.8|l now simplify to 

R'l = R l [-(eAl + £A*) 2 + (eA* + £A*) 2 ]+2\ 1 RAR 2 + R 2 -u 2 ) + \ 2 R l R 2 , (4.10) 
A\ x " = 2eR 2 {eA\ x + eA*' x ) + 2eR 2 (e A*' x +eA\ x ) (4.11) 

The final choice is whether to choose the gauge fields to be timclike (A\ 0,Af — 0), lightlike (A\ — Af ^ 0), or 
spacelike (A\ = 0, Af 7^ 0). As is obvious from Eq. (|4.1U|I . the timelike ansatz gives rise to asymptotic oscillatory 
behaviour for the R's and is thus unacceptable. We first discuss the spacelike case. 

As a boundary condition, we require that i?i i2 tend to vacua 1 and 2 on opposite sides of the wall: 

i?i(+oo) = u, R 2 {+oo) = and R^-oo) = 0, R 2 (-oo) = u. (4.12) 

The appropriate boundary conditions for the gauge fields are then determined by examining Eqs. i|4~TT|) at the 
asymptotic points. As w — > +00, we find that 



I e ! «' t , (4.13) 



^1 

A 2 " I \ ee e 2 I \ A 2 



where Aj = Af . Inputting trial asymptotic solutions of the form 

A h2 =a lt2 e- KW , (4.14) 

we obtain an eigenvalue equation yielding 



K = V2u \J e 2 + e 2 and k = 0. (4-15) 

The zero eigenvalue means that one linear combination of A\ and A 2 is always unsuppressed, and the vanishing second 
derivative tells us it grows linearly with w (provided its amplitude is nonzero). By the symmetry of the problem, a 
similar situation obtains as w — > — 00. The exponentially suppressed linear combination embodies the Meissner effect. 
Numerical solutions are depicted in Fig. El 
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FIG. 5: Static solutions to the two scalar and two gauge fields in the U(1)®U(1) model. In all cases Ai = A2 = e = 1, 
and the fields are all plotted in units of u. In the top left plot, e = as in the Rozowsky et al. model and the gauge fields 
show Meissner-like suppression under their respective scalar field. In the other three plots, e = \ and we show cases with 
various values for the free gauge field boundary conditions. In the two right plots, reflection symmetric boundary conditions are 
imposed for the pair of linear combinations eAf + eAj . In the bottom left plot, the gauge boundary conditions are asymmetric 
and their coupling to the kink shifts it to the right. 



The lightlike case is qualitatively similar to the spacelike case. But it is amusing to point out that for the parameter 
choice A2 = 4Ai an analytic solution exists. For the lightlike ansatz, the gauge fields cancel out of Eg. 14.11)1 and with 
the stated parameter choice and e = 1, e = 0, the resulting equations yield the solutions 

u r — 

R h2 (w) = - |l ± tanh(VAi «w)J . (4.16) 

Inputting these results into Eq. 14.111 the analytic solutions obeying the correct boundary conditions for the gauge 
field functions A\ = Af = Aa are 



Ai, 2 (w) = const, x a F x ( , I+y|WI . 1 + ^2 ; (1 + e* 2 ^™)- 1 ) (l T tanh^u w)) 75 , (4.17) 

where 2F1 is a hypergeometric function. 
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B. Warped spacetime 

We now generalise the above model through the inclusion of gravity 13] . The action is 

S = -[[%R+ g MN t MN + \g MP g NQ f M NPQ + v]^—g d 5 x, (4.18) 

with 

t MN = ^{D M $i)* D N $i, (4.19) 

i 

D M = d M -^iQiA iM , (4.20) 

i 

hiNPQ = FjMNFjPQ, (4-21) 

i 

and the charges (Qi, Q2) of the scalar fields are (e, e) for $1 and (e, e) for $2- 
The Higgs and gauge field equations of motion are 

F)V 

Dm(V^9 9 MN D n ^)-V^9^-; = (4.22) 

Qm (V=9 9 MP 9 NQ FiP Q ) +V=g~g NP 2 (eIm(^Dp*<) + elm (tyDpfy)) = 0, (4.23) 

while the Einstein equations, kGmn — —Tmn, feature 

Tmn — 2tMN + g PQ Jmpnq + gMN&U, (4.24) 
Cm = -9 MN t M N - \9 MP 9 NQ fMNPQ V. (4.25) 

The metric ansatz is generalised to 

ds 2 = -e-^dt 2 + e h ^dx 2 + e^ w) (dy 2 + dz 2 ) + dw 2 , (4.26) 

in order to have enough degrees of freedom to be able to look at sufficiently general gauge field configurations. We 
look for solutions of the form 

A m = (Ai{w), Bi(w), 0,0, ZiM), (4.27) 
$i = R l (w)e ia ^ w \ (4.28) 

To simplify the algebra it is useful to define the following: 

(4.29) 
(4.30) 
(4.31) 
(4.32) 
(4.33) 
(4.34) 

In terms of these, the equations of motion for the gauge field components are 

A'l + 7 ' A' i -2eR?A i -2eR 2 j A j = 0, (4.35) 

B" + y' B[ - 2eR 2 i B i ~ 2eR?Bj = 0, (4.36) 
(ei? 2 + ei? 2 ) Z t = 0. (4.37) 



T = 


W 


+ h + 


2i), 


y = 


K- 


f + h 


+ 2j) , 


7 = 


W 


-h + 


2i), 




eAi 


f eAj, 






eB i 


f eBj, 




Zi = 


eZi - 
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Equation i|4.37[l implies Zi = which we use to simplify the rest of the equations of motion. For the scalar fields we 
have 

8V 

R'l + T'R\ + e~ f A 2 ,R t - eT h B?Ri - e~ ia * ^ = 0, (4.38) 

and for the three metric functions we have 

f" + if 2 + \fti + iff \h'f \f - ge"/* - + ±n = o, (4.39) 

h" + ifh' iff + \h a + Ih'f - \f + |c^$ + §a e -"* + A fi = o, (4.40) 

f ~ y'ti + hf'f + ih'f + If + J-e-'* - ^c- h * + §-n = 0, (4.41) 



where 



$ = AlRl+AjRl + ^Af + ±A'i, (4.42) 



2' 

* = B\R\ +BlRl + \B? + \B'i, (4.43) 
Q = Rf + Rf + V. (4.44) 

We performed an extensive numerical search for solutions to Eqs. (|4.35l |4~3U1 14.38I4.41|I . subject to the boundary 
conditions that the Higgs fields tend to VEVs, and that the metric functions tend to Randall-Sundrum form. For 
dcfiniteness, we used the same Higgs potential as in the non-gauged model, namely as given by Eqs. I|3.20|l . (|3.22|) 
and (|3.25|l . The result of this investigation was the perhaps disappointing one that all solutions found required the 
gauge fields to vanish. Since the investigation was numerical, we cannot be certain that no non-trivial solutions exist, 
although we strongly suspect this to be the case. 



V. CONCLUSION 

We have constructed a U(1)<8>U(1) model with two complex scalar fields interchanged by a discrete symmetry 
coupled to 4 + 1 dimensional gravity. It has a solution featuring clash-of-symmetries-style Higgs kink configurations 
in a Randall-Sundrum-like spacetime. Gravity is localised to the dynamically generated (smooth) brane, and the 
Randall-Sundrum limit of the solution can be defined. For the chosen Higgs kink and metric configurations, the 
Higgs potential had to be of a certain sextic form whose properties we studied in some depth. The symmetry breaking 
pattern varies as a function of the extra dimension coordinate w, and displays the clash of symmetries phenomenon: at 
all points \w\ < oo both U(l)'s are broken, with alternate U(l)'s restored asw-> +oo and w — > — oo. The spontaneous 
breaking of the discrete interchange symmetry guarantees topological stability for the Higgs-gravity-induced brane. 
We also examined the gauged version of the theory, but we found that to have nonzero gauge fields the spacetime has 
to be flat. 

This work sets the stage for incorporating gravity into more complicated models displaying the clash of symmetries. 
The eventual aim is to construct a realistic brane-world model that uses the clash of symmetries to understand 
spontaneous symmetry breaking in a more satisfactory fashion than obtained with the default mechanism using 
homogeneous Higgs vacuum expectation values. 
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